Local density of states at zigzag edge of carbon nanotubes and graphene 
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The electron-phonon matrix element for edge states of carbon nanotubes and graphene at zigzag 
edges is calculated for obtaining renormalized energy dispersion of the edge states. Self-energy 
correction by electron-phonon interaction contributes to the energy dispersion of edge states whose 
energy bandwidth is similar to phonon energy. Since the energy-uncertainty of the edge state 
is larger than temperature, we conclude that the single-particle picture of the edge state in not 
appropriate when the electron-phonon interaction is taken into account. The longitudinal acoustic 
phonon mode contributes to the matrix element through the on-site deformation potential because 
the wavefunction of the edge state has an amplitude only on one of the two sublattices. The on- 
site deformation potentials for the longitudinal and in-plane tangential optical phonon modes are 
enhanced at the boundary. The results of local density of states are compared with the recent 
experimental data of scanning tunneling spectroscopy. 
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I. INTRODUCTION 

The local electronic property near the edge of graphene 
depends on the lattice structure of the edge. For example, 
the zigzag edge induces edge states which are 7r-electron 
states localized near the edge while the armchair edge 
does noti The edge states which have a flat energy dis- 
persion, show a peak in local density of states (LDOS) 
near the Fermi energy. The peak structure in LDOS has 
been observed at the edge of graphene by scanning tun- 
neling microscopy (STM) and spectroscopy (STS)^^^ 
The LDOS peak is a direct evidence of the edge states be- 
cause the peak is not observed at an armchair edge and 
the height of the peak decreases with increasing a dis- 
tance of the STS tip on graphene plane from the zigzag 
edge. The data on LDOS are useful to understand the 
energy and life-time of an electron in the edge states. 
The life-time of the electron is determined by electron- 
phonon (el-ph) interaction and the el-ph interaction is 
important for almost flat energy dispersion when the De- 
bye energy (Hcud ~ 0.2 eV) is comparable to the energy 
bandwidth. Thus the el-ph interaction for edge states af- 
fects STS spectra and is essential for an analysis of STS. 
In this paper, we consider self-energy correction for edge 
states induced by the el-ph interaction, and compare the 
theoretical results of LDOS with experimental data^^ 

A complete flat energy dispersion relation of the edge 
states is widely recognized by the theory^ However, 
STM/STS2i2iii§i and angle-resolved photo-emission spec- 
troscopy (ARPES) 6 show that the edge states have a 
small but finite energy dispersion. Using STM/STS at 
graphene edge, Niimi et al£& and Kobayashi et al£& 
independently observed a peak in the LDOS below the 
Fermi energy by 20 ~ 30 meV. The peak position rela- 
tive to the Fermi point (Ep — 0) shows that the edge 



states have a finite bandwidth. Using ARPES, Sugawara 
et air- observed the Fermi surface of Kish graphite and 
found a weakly dispersive energy band near the Fermi 
energy. In the previous paper, we pointed out that next 
nearest- neighbor (nnn) tight-binding Hamiltonian, Tl nnn , 
is essential for the bandwidth.- As shown in Sec. HH 
Hnnn lowers the energy dispersion of the edge states as 
E(k) = 7„(2 cos ka + 1) (2ir/3 < ka < 4ir/3) where j n 
and a is the hopping integral between nnn sites (7„ ~ 0.3 
eV) and a lattice constant of graphite (a = 2.46 A), 
and the value of 7„ is calculated on the basis of density- 
functional theory by Porezag et al& Since ka = ir state 
is located at the bottom of the band and ka — ► 27r/3 (or 
47r/3) is located at the top of the band, the bandwidth 
of the edge states, W, is given by W — 7„. However, 
the observed energy bandwidth (20 ~ 30 meV) is much 
smaller than 7„. The reason why observed bandwidth is 
smaller than j n is that the self-energy correction 
renormalizes E(k) as E(k) + Re(£(fc)). Because the el- 
ph interaction makes the effective mass of the edge states 
large, W generally decreases by taking account of el-ph 
interaction. It is also noted that T(k) = — 2Im(£(fc)) rep- 
resents the energy uncertainty of the edge state. Since 
the Fermi-Dirac distribution function has a width of k^T 
around the Fermi energy (Ep), if T(k) > k^T, then it is 
not appropriate to treat the edge state by a single-particle 
picture. 

The el-ph interaction is calculated by the matrix ele- 
ment of deformation potential. When the wavefunction 
is expanded by tight-binding orbitals, the matrix clement 
consists of on-site and off-site atomic deformation poten- 
tials £ The el-ph matrix element of a given wavefunction 
is given by the sum of atomic deformation potentials over 
all carbon sites on which the electronic wavefunction has 
an amplitude. The el-ph interaction for edge states shows 
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a different behavior from that for extended states. The 
unit cell of graphene consists of two sublattices; A and 
B. The extended states have a finite densities on both 
sublattices while the edge states have density only on 
one sublatticei Suzuura and Ando pointed out for the 
extended states that the on-site atomic deformation po- 
tentials at A-atom and B-atom cancel with each other in 
the matrix element for a backward scattering process be- 
cause of a phase difference of the wavefunctions for two 
sublattices Thus only the off-site atomic deformation 
potential, that is generally weaker than the on-site one, 
contributes to the backward scattering. This is consis- 
tent with the fact that a metallic carbon nanotube (CNT) 
shows the quantum conductance and ballistic character 
at a low temperature ! 11 ' 12 ? 13 However, the cancellation 
of on-site deformation potential does not occur for edge 
states, since the wavefunction has an amplitude only on 
a sublattice for the edge state. Thus we can expect a 
relatively strong el-ph interaction for the edge states and 
a large self-energy correction to the edge states. 

The el-ph interaction for the edge states is relevant to 
many observations in experiment of CNTs. Superconduc- 
tivity in CNTs is an important example. Takesue et al^ 
observed a drop of resistivity in multi-walled CNTs and 
pointed out that the connection of multi- walled CNTs to 
Au-electrode is sensitive for the occurrence of the resis- 
tivity drop. Since the edge states enhance LDOS near 
the ends of a CNT, the el-ph interaction should be sensi- 
tive to the properties at the interface between the CNT 
and an electrode. Furthermore, we propose that the large 
LDOS and a strong el-ph interaction favor superconduct- 
ing instability for the edge states.— The self-energy cor- 
rection is important for an estimation of the supercon- 
ducting transition temperature. Thus a quantitative dis- 
cussion of the el-ph interaction for edge states will play a 
decisive role for a future work on STS and superconduc- 
tivity of graphene based materials. 

This paper is organized as follows. In Sec.|TT]we give a 
method for calculating el-ph interaction for edge states. 
In Sec. IIII1 we show phonon mode and localization length 
dependence of el-ph interaction. In Sec. IIV1 we calculate 
self-energy correction to edge states, by which we obtain 
a renormalized energy dispersion relation and LDOS. The 
theoretical results are compared with the experiment. 
Discussions and summary are given in Sec. |Vj 



II. ELECTRON-PHONON INTERACTION FOR 
EDGE STATE 

A. Edge States 

The edge states of graphene are 7r-electronic states lo- 
calized near the zigzag edge (Fig. QTa)). 1 Since a CNT is 
a rolled-up sheet of graphene, the edge states may exist 
near the open boundaries of (n, 0) zigzag CNT regardless 
of the value of n (see Fig. []Jb)). The energy eigen equa- 
tion of the nearest-neighbor (nn) tight-binding Hamil- 



tonian, W n n|*fc) = E{k)YS>k), gives a flat energy band, 
E(k) — 0, for the edge states where k is the wavevector 
along the edge. As shown in Fig.QJc), the edge state ex- 
ists when 2n/3 < ka < 4ir/3 and 7i nrm lowers the energy 
dispersion of the edge states as E(k) — W(2 cos ka + 1) 
(2ir/3 < ka < 4ir/3)Z The edge state behaves as a plane 
wave around the axis (R c in Fig. [TJb)), while wavefunc- 
tion is localized in the direction of nanotube axis i? 4 . 
The localization length in the direction perpendicular to 
the edge is fc-dependent, £(fc) = -|T|/2 ln(2| cos(fca/2)|) 
(see Fig. (Ud)). 16 Though £(fc) becomes infinite when 
ka — > 27r/3 or ka — > 47r/3 for a graphene, we can show 
that £(fc) < dt/2 for a CNT where d t = na/ir is the 
diameter of the zigzag CNT. To explain this, let us con- 
sider a metallic zigzag CNT, (3m, 0). Then k for the 
edge states are discrete due to the periodic boundary 
condition, which is given by k(i) — 2n/3a + 27ri/3ma 
(i = 1, . . . , to — 1). We get the largest localization length 
£(*(!)) =£(fc(m-l))« dt/2. 




FIG. 1: (a) The unit vectors of graphene are denoted by 
ai and &2. ai is the unit vector around the tube axis and 
T (= 2a2 — ai) is the translational vector along the zigzag 
tube axis. We define 21 = |T| where I = 0.21 nm and a = 
ai|. The density of an edge state has a value only on A- 
atoms, which is represented by solid circles. The radius of 
a solid circle is proportional to the density, which shows the 
localization, (b) The edge states exist near the zigzag edge 
(7?* = 0) of the (n, 0) CNT with open boundary. The diameter 
and length of the CNT is denoted by d t and L, respectively, 
(c) Tinn has a flat energy band of the edge states between the 
K and K' points at the Fermi energy (E-p — 0). Tinnn causes 
a bandwidth W — 7„. (d) We plot the localization length, 
£(k)/l = -l/ln(2| cos(fco/2)|) for 2tt/3 < ka < 4tt/3. 

The wavefunction of the edge state is written as 7 

I**) = V^E eX P \* kR u,A- fir) l#Ru,A)>, (1) 
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where |</>(R U]S )) is 2p z orbital of a carbon atom and N k is 
a normalization constant. The summation on u is taken 
for all unit cells of graphene or a CNT. We take the co- 
ordinate R = (i? c , i?') on the cylindrical surface of a 
zigzag CNT, in which R c and R* are coordinates around 
and along the tube axis, respectively (Fig.QJb)). The po- 
sition of a carbon atom is denoted by R p = s , -R„ s ) 
where p = (u, s) represents the s-th sublattice (s = A, B) 
in the u-th hexagonal unit cell. As is taken for the 
zigzag edge site iZ* A =0, the edge state has amplitudes 
only on A-atoms (s = A). Equation |T]) is correct for 
27r/3 < ka < tt. For tt < ka < 47r/3, a phase factor 
exp{i7ri?^ A /£} should be multiplied to |</>(R UiJ 4)) J£ 



B. Electron-Phonon Interaction 

The el-ph interaction for graphene is formulated by 
Jiang tt aZ.— , which will be applied to the edge states. 
The el-ph interaction for the edge states is written as 



rtint — 
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:]r^c^,(q)(& q ,, + &L q >i, Cfc , (2) 



k,k' qt,v 



where N u is the number of graphite unit cells, c k is the 
annihilation operator of the edge state, and b^ v is the 
annihilation operator of the v-ih phonon mode. There 
are six phonon modes: out-of-plane tangential acous- 
tic/optical mode (oTA/oTO), in-plane tangential acous- 
tic/optical mode (iTA/iTO), and longitudinal acous- 
tic/optical mode (LA/LO).™ a£ fe , (q) is the el-ph in- 
teraction connecting two edge states k and k' by t^-th 
phonon mode with momentum q. Due to the momen- 
tum conservation along the edge, k' = k + q, while the 
wavevector perpendicular to the edge qt is needed to sum 
over the Brillouin zone. a kk , (q) is given by a kk , (q) = 
A w (q)M£ fc ,(q)/\/2 where A v (q) = y/K/m^oAc^ is the 
amplitude of phonon (7ki;„(q) is the energy of the t^-th 
phonon with the momentum q) and M kk , (q) is the el-ph 
matrix element, 



M&,(q) = -^<^|V«(Rp)|* fc ) • t/(R p )e q (s) e ^ R P . 



(3) 



Here i>(R p ) is the Kohn-Sham potential of a neutral pseu- 
doatom calculated on the basis of density-functional the- 
ory by Porezag et al& for a carbon atom at R p , e q (s) is 
phonon eigenvector at an s-atom normalized in the unit 
cell as J2 S =A b e q( s )* ' e q ( s ) — 1j an d U(R P ) is a rota- 
tional operator for e q (s) from an s-th atom at origin to 
a s-th atom at R p . To obtain w„(q) and e q (s), we use 
the force constant parameters calculated by Dubay and 
Kressei^ for the dynamical matrix^ 



Putting Eq. (TJ into Eq. Q, we obtain 

N k >N k 



M£ fc ,(q) = - 
X eXp \ ~ ik ' R n',A + ik K,, 

u' ,U \ 

x m(R„',A, R«,a; q, v), 



(4) 



where m(R M ' jJ 4, R Mj a; q, v) is the atomic deformation po- 
tential, 9 defined by 



m(R u /,A,R u ,A;q J v) 

^(R u%a )|V^(R p )|^(R u ,a)) • U(R p )e^( s )e^ R -. 

(5) 



There are two types of m(R u ' iS /, R U]S ; q, v). The first 
type is the case of (it', s') = (it, s) which is referred to as 
the on-site atomic deformation potential. 9 The other one 
is (it', s') 7^ (it, s) which is the off-site atomic deforma- 
tion potential. The on-site (off-site) atomic deformation 
potential represents a scattering process of an electron 
from Ru, s to the same (different) site. 

The off-site deformation potential matrix element for 
the next-nearest A atoms, |(0(R u / :J 4)|Vu(R p )|0(R tt ^)) • 
n\, is negligible (for any R p ) in Eq. (f5]) because it de- 
cay quickly as a function of |Rv ti 4 — R«',a| where h 
is a unit vector along the two carbon atoms. Density- 
functional theory gives that the off-site deformation po- 
tential matrix element for nearest-neighbor interaction 
is |(^(R u ,b)|Vt;(R u ,b)^(R u ,a)) • n| ~ 3eV/A. 8 - 9 As we 
pointed out above, the wavefunction of the edge state 
has an amplitude only of one sublattice and thus this 
nearest- neighbor term does not contribute to a kk , (q). 
Thus the off-site atomic deformation potential does not 
contribute to a kk , (q) for the edge states. For the on- 
site atomic deformation potential, we need to consider 
several carbon atoms which are located near R^.^ for 
the center of deformation potential R p in Eq. (JSJ). The 
value of |(^(R UiA )|Vu(Rp)|^(R u , A )) • n| is not negligi- 
ble if |Rp — Ru,a| < 3 A. Density- functional theory gives 
that the largest contribution from nearest-neighbor site is 
|(0(R„,a)|V W (R^ s )|0(R u , a )) ■ n| ~ 8 eV/A. 8 Thus on- 
site deformation potential is more important than the 
off-site deformation potential^ 

Now we can write m(R u ^A, H U ,A] q, v) as 

m(R u ' ij4 ,R Uij4 ;q, v) = <5„'„m on ( J R^ A ;q,j/)e lq - R "- 4 , (6) 

where 77i on (i?^ A ; q, v) is on-site deformation potential 
from all possible R p , defined as 

TO on(H* U!A ;q, v) = 

^<0(R m , a )|Vu(R p )|^(R„,a)) • C/(R P )e q (s)e^( R -- R -^. 
p 

(7) 
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Putting Eq. © into Eq. (gj), we get 

M fc V(q) = -N k ,N k 6 k ,, k+q 

(8) 

For ka < n < k'a or k'a < tt < ka, we must multiply 
expliirR^ ^/i} to m on (i?* A ; q, i/) in Eq. ((HJ. 

It is noted that m on (R t u A ; q, f) is defined to be inde- 
pendent of i?£ A . Instead, R^ A appears in the phase 
of to(R u . j 4, R n .yi; q, f) in Eq. ((6]). The phase gives the 
momentum conservation around the tube axis, S k ',k+qj 
after the summation about R c u A is made in Eq. Q . On 
the other hand, m on (i?^ A ;q,i>) depends on R l u A because 

J2 P m Eq. is restricted for i?p > and m on (-Ru,A; q> v ) 
does not have a translational symmetry near the edge. 
Hereafter, we refer m on (ii* A ; q, as boundary defor- 
mation potential. 

Won a ; q, v) depends on v strongly. The iTO and 
LO modes whose eigenvectors are pointing in the direc- 
tion perpendicular to the edge give a large contribution 
to m on (R t u A ;q, v). We generally expect that the bound- 
ary deformation potential appears only near the edge. In 
fact, |(0(R U) a)|Vu(R p )|^(R„ 1 a)) • n\ is finite only when 
|Rp— Ru^l < 3A. Thus m on (i?^ j4 ; q, v) depends on R* u A 
only when R l u A < 3A. The boundary deformation poten- 
tial contributes to the matrix element increasingly with 
decreasing or which can be seen in Eq. {§]). 

In Eq. I|7|). we must consider relative atomic motion of 
R U ,A to Rp. This can be done by replacing U(R P ) with 
U(Rp) — E where E is the unit matrix. In Fig. Efa) we 
show phonon eigenvector of the oTA mode pointing per- 
pendicular to the nanotube surface. If we consider the 
phonon eigenvector in a flat surface or graphene, the rel- 
ative displacement of the two atoms becomes zero when 
q = 0. However, it is not the case for a cylindrical sur- 
face. As shown in Fig. HJb), the relative eigenvector, 
([/(R p ) — £')e° TA (s), remains finite. The relative vec- 
tor for the oTA mode changes the area on a cylindrical 
surface and is similar to the LA mode except for a l/d t 
reduction of the amplitude, by which m on (i?^ A ; q, oTA) 
is proportional to l/dt- On the other hand, the relative 
displacement becomes zero for the in-plane modes such 
as iTA and LA, when q — > 0. Thus, the el-ph interaction 
by the iTA and LA modes vanish for q = while the 
oTA mode provides a finite m n(Rt A \ q, oTA) even for 
q = 0. Hereafter, we simply use U(R P ) for U(R p ) — E. 



III. CALCULATED RESULTS 

In this section, we plot \M% k , (q)| for several values of 
k and k', and examine the dependence of |M£ fe ,(q)| on 
phonon mode v in Sec. IIII II and on nanotube diameter 
dt in Sec. InOI 



(b) (U(B P )-E)e^ A (A) 

e% TA (A) U(H p )e^ TA (A) relativevector 




FIG. 2: (a) For the oTA mode with q = 0, the eigenvectors of 
two atoms at R u ,a and R p are pointing perpendicular to the 
CNT surface, (b) In order to calculate the atomic deformation 
potential, one needs to consider the relative displacement of 
two carbon atoms. 



1. v dependence of \M kk ,(q)\ 

In the calculation, we consider the magnitude of 
M£ fc ,(q) for the (60,0) CNT (d t » 5 nm). In Fig. 
we plot \M% k ,(q)\ as a function of q t for (a) (k, k') = 
(77r/10,87r/10) and (b) (26tt/30, 29tt/30). The corre- 
sponding phonon eigenvector is the same (q = 7r/10). 
Thus the difference between the two cases shows the de- 
pendence of |M£ fc ,(q)| on and £(*;'). (7tt/10, 8tt/10) 
is chosen as an example that £(fc) (= 22A) and 
(= 4.4A) are longer than the carbon-carbon bond length 
(a cc = 1.4A), while (267r/30, 297r/30) is chosen as an ex- 
ample that £(fc) (= 2.4A) and (= 0.9A) are compa- 
rable to a cc . 

As shown in Fig.^a), |M™(q)| is 4 eV/A at q t = 0. 
|M^,°(q)| decreases with increasing q t , while |Mjy£?(q)| 
increases with increasing q t . This behavior of the iTO 
and LO modes relates to the boundary deformation po- 
tential. The eigenvector of the iTO (LO) mode is point- 
ing along the tube axis when g*/|T| < <?/|ai| (<?t/|T| > 
q/|ai|) and then produces a large boundary deformation 
potential. On the other hand, e° TO (s) is perpendicu- 
lar to the CNT axis and the oTO mode does not con- 
tribute to a boundary deformation potential. The value 
of |M£j,°(q)| is considerably smaller than |M£ fc ,(q)| for 
the iTO and LO modes. Thus the contribution of the 
oTO mode to the el-ph interaction can be neglected for 
rf t «5nm CNT. In Fig.^b), |M^,°(q)| and |Af fc L fc 9(q)| 
reach 11 eV/A, which indicates that |M£ fc/ (q)| for the 
iTO and LO modes increase significantly with decreas- 
ing and £(fc'). The boundary deformation potential 
becomes more effective with decreasing the localization 
length. 

The matrix element for iTA mode is smaller than the 
other acoustic modes for a wide range of qt as shown 
in Fig. [3ja). Though \M^, (q)| can be comparable to 
\MuB (q)| as shown in Fig.[3jb), the contribution of the 
iTA mode to the el-ph interaction is negligible to that for 
oTA because the amplitude of the iTA mode is smaller 
than that of the oTA mode; A iTA (q) < A oTA (q). On 
the other hand, the oTA and LA modes are important 
for lower temperature in the el-ph interaction. The oTA 
mode changes the volume of a CNT and gives an on- 



(a) 



(k,k') = (77r/10,87r/10) 
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(b) 



(k,k') = (267r/30,297r/30) 




FIG. 3: (Color online) |M£ fc ,(q)| for the (60,0) zigzag CNT: (a) {k,k') = (7tt/10, 8tt/10) and (b) (k,k') = (26tt/30, 29tt/30) 
are plotted as a function of qt where q = tv/10. Three solid curves represent acoustic phonon modes: oTA(green), iTA(blue) 
and LA (red), and three dashed curves are optical modes: oTO (green), iTO(blue) and LO (red). The vertical dashed lines 
represent q t = V3q (y/S = |T|/|ai|). 



site deformation potential as shown in Fig. [2l Moreover, 
the energy of the oTA mode is smallest among acoustic 
modes and thus A oTA (q) can be larger than A LA (q). The 
LA mode is a area-changing mode and produces a large 
on-site deformation potential, and contributes to the ma- 
trix element most significantly among acoustic modes. 



2. dt dependence of \M£ k , (q)\ 

It is shown from Eq. © that M% k , (q) for two CNTs 
with different diameters (dt) are the same for the same 
values of k and k'. In general, we can find the similar 
values of k (or k 1 ) for different dt- However, it is not the 
case for edge states near the K or K' point. Then M kk , (q) 
depends on d t for such edge states. For instance, M kk , (q) 
depends on d t for the largest £(k) « dt/2. In Fig. 0J we 
plot |M£ fe ,(q)| for (41tt/60, 47tt/60) of the (120,0) CNT 
(d t ~ 10 run). \M kk , (q)| has the similar functional shape 
as shown in Fig. [3^a) while the values become smaller. 

The reduction of \M kk , (q)| can be explained as follows. 
If we neglect dt dependence of the boundary deformation 
potential, the summation on R u A in Eq. ([5]) can be per- 
formed analytically. M kk , (q t = 0) is then proportional to 
the factor Rkk'{dt)', 



Rkk< {dt 



(9) 



where Rkk'(dt) is a function of d t because £ depends on 
dt- Since Rkk'(5nm) = 0.37 and i?fcfc'(10nm) = 0.30, we 
have R kk ' (lOnm) / R kk i (5nm) w 0.81. This ratio repro- 
duces |M fc L A( q )| 10nm /|M fc L A(q)| 5nm = 1.3/1.5 « 0.86 at 
qt = 0. 



(k,k') = (4l7r/60,477r/60) 




FIG. 4: (Color online) |M£ fc ,(q)| for (120,0) zigzag CNT. The 
phonon eigenvector is the same as Fig. [3] The vertical dashed 
line represents qt = V%q. 



IV. LOCAL DENSITY OF STATES 

Now we calculate a renormalized energy dispersion re- 
lation and LDOS using self-energy £(fc, iui n ) which is de- 
fined by 



EE 



Nu ^ ^ ( w « ~ w ™) 2 + (^(q)) 2 



iuj m - (E(k + q) - E F ) - £(fc + q, iui m ) ' 

(10) 
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where T is temperature and ui n — irkBT(2n + 1) is the 
Matsubara frequency (n is integer). We put the cut-off 
Matsubara frequency, \ui n \ < hup, ~ 0.2 eV. By means of 
Pade approximation, 19 Y*(k,iujn) is changed to that on 
the real-axis of uj: E(fc,cj). Then we find a solution e 
satisfying e = {E(k) — Ep) + Re (S(/c, e)) as a function of 
k, Ep and T. The obtained CE F ,T{k) is a renormalized 
energy dispersion relation. We adopt T = 50 K in the 
following calculations. The T-dependence of eE F ,T(k) is 
negligible for a wide range of T, for instance, T — 77 K 
gives almost the same result. Calculations at low temper- 
ature have some numerical advantage since the number 
of u n increases. The LDOS curve is defined as a function 
of bias voltage (V) between graphene and the STS tip, 
and the distance (iZ*) from the zigzag edge sites and the 
STS tip by 

(ii) 

where r BF (fc) = -2Im(E(jfc, e EF (k))) is the width of STS 
spectrum. It is noted that the value of Ep is given in 
such a way that the number of valence (edge) states is 
conserved when we calculate E(fc,iw„) self-consistently, 

that IS, Y,E(k)<E F 1 = Y,e EF (k)<E F 1- 

A. Numerical Result 

In Fig. Eta), we plot E(k) (W = 0.3 eV) without 
self-energy correction as the dashed curve and renormal- 
ized energy dispersion e£ F (fc) for Ep — —0.054, —0.122, 
—0.250 eV as the blue, red and green curves, respec- 
tively. The reason why we consider different values 
is that in the experiment charge transfer from STS tip 
or substrate may modify the Ep values and that we need 
to investigate Ep dependence. The parallel dashed lines 
denote the Fermi level for these Ep values. For each 
value of Ep, Im(E(fc, es F (fc))) is plotted to show the k- 
dependence. By denoting the renormalized bandwidth 
as W, we obtain W w 0.2 eV for Ep = -0.054 eV 
and —0.122 eV. The corresponding mass enhancement 
parameter A is about 0.5 where we defined A = W/W' — l. 
When Ep = -0.250 eV, W is about 0.27 eV and A « 0.1. 
We observe that the value of TE F {k = n) is about 0.13 
eV when Ep = —0.054 eV. The calculated e EF (k) and 
Im(E(fc, €E F (k))) strongly depend on the value of Ep. 

The values of W and r# F (fc) relate to the peak posi- 
tion and width of D(V, R l ). We plot D{V, i?*) for Ep = 
-0.054 eV at i?* = 0.42, 0.84, 1.26 nm in Fig. [3b). The 
LDOS curves have several sharp spikes at V = €e f (k) due 
to relatively small value of r^ F (fc) compared with the fi- 
nite level spacing of the edge states. Here we take (120, 0) 
CNT (dt ~ 10 nm) for calculating the self-energy. For a 
graphene (dt — > oo), the level spacing becomes zero and 
the spike structure will disappear. The calculated LDOS 
structure for Ep = —0.054 eV has a peak near V « —60 



meV which is close to the observed LDOS in which a peak 
is located at V = —30 ~ —20 meV.->&^ However, the 
width of the peak is about 0.7 eV, which is much larger 
than the experiment 2,4 (0.05 ~ 0.1 eV). The peak posi- 
tion and the width are improved to fit the experimental 
data when W — 0.1 eV as shown in Fig. (5^c) and (d). 
In this case, W = 0.02 and 0.04 for Ep = -0.002 and 
—0.27 eV, respectively. The calculated LDOS structure 
for Ep = —0.027 eV has a peak near V « —20 meV with 
the width of ss 0.05 eV. The small value of W does not 
show any spike structure due to the discrete k values. 
Thus overall feature seems to be better for W = 0.1 eV 
than W = 0.3 eV. It should be noted that W = 0.1 eV 
is not always consistent to the experiment. When we as- 
sume all the edge states are below the Fermi energy, we 
see that Ep appears above E = eV which gives a peak 
very close or above the V = eV, while the experiment 

shows V = -30 20 meV. 

It is worth mentioning that Niimi et al.—^ reported that 
the tunnel current was unstable when the STS tip was 
located very close to the zigzag edge. Then the electron 
injected from the STS tip has a large transition ampli- 
tude to edge states having £(fc) < 21 (0.42 nm), that is, 
fc-states which are around k = tt state (see Fig. Hfd)). 
As shown in Fig.OJc), the magnitude of TE F {k) is much 
larger than k^T (w 0.0045eV) for most value of k and 
r.E F (fc) for states around k = it state reaches about 0.02 
eV and yield strong fluctuation. It indicates that the 
tunnel current is unstable. We calculate D(V, 0) and find 
that the height and width of the peak are both signifi- 
cantly larger (more than 10 times larger) than D(V, i?*) 
for R l = 0.42 nm. As we noted in Sec. M iTO and LO 
modes give a large matrix element through the bound- 
ary deformation potential. The boundary deformation 
potential may be relevant to the unstable tunnel current. 
Since the injected electron from a STS tip is localized 
near the edge, we expect that the tunnel current would 
be strongly affected by the boundary deformation poten- 
tial. 



V. DISCUSSION AND SUMMARY 

We showed that el-ph interaction for edge states con- 
sists only of the on-site atomic deformation potential. As 
a result, LA mode contributes to scattering most effec- 
tively and the on-site deformation potential is enhanced 
at the edge for LO and iTO modes. It is to be noted 
that the on-site atomic deformation potential does not 
contribute to backward scattering of extended states and 
the off-site atomic deformation potential gives rise to re- 
sistivity^ Because the on-site atomic deformation po- 
tential is larger than the off-site atomic deformation po- 
tential, the edge states exhibit the strong el-ph coupling 
character that the graphite system originally possesses. 

The original energy bandwidth of the edge states, W, 
is consistent with the observed position of LDOS peak 
when W w 0.1 eV, which is the same order of 7„ w 0.3 
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FIG. 5: (Color online) Calculated energy dispersion e Ep (fe) and LDOS D(V, R*) for (a,b) W = 0.3 eV and (c,d) W = 0.1 eV. In 
(a), E(k) is plotted as the dashed curve, and eE F {k) is plotted for Ef = —0.054, —0.122, —0.250 eV. We plot Im(S(fc, eE F (k))) to 
show the fc-dependence. In (b), D{V, R l ) for Ef = —0.054 eV and its smoothed curve (dashed curve) are plotted for R l = 0.42, 
0.84 and 1.26 nm. In (c), EE ¥ {k) is plotted for Ef = -0.002 eV and -0.027 eV. In (d), we plot D(V, fl*) for Ef = -0.027 eV. 



eV but not the same value. It is noted that the case of 
W = 0.3 eV does not include the overlapping integral 
(s parameter—) which increases (decreases) the conduc- 
tion (valence) bandwidth. To examine the effect of s 
parameter on the bandwidth of the edge states, we per- 
formed the energy band structure calculation in an ex- 
tended tight-binding framework^ and obtained W « 0.2 
eV. We expect that W is externally modified by attach- 
ing a functional group or contacting an electrode to the 
edge sites. In addition to the el-ph interaction, electron- 
electron (el-el) interaction causes self-energy which may 
account for the difference. The cl-cl interaction con- 
tributes to Yi(k,iuj n ) additively when there is no cross 
term of el-ph and el-el interactions. Since many papers 
report the importance of el-el interaction, especially the 
contribution to the imaginary part of self-energy should 
be important and thus the spike structure of discrete k 
does not appear in the real case. Then W decreases 
and the width of LDOS increases. The effect of el-el in- 
teraction warrants future work concerning the effects of 



dynamical details on LDOS curve. A detail experimental 
data of STM/STS for zigzag edge of CNTs and graphene 
may be useful for a qualitative estimation of the strength 
of the el-el interaction. When the values of W and Tilo-o 
are comparable, the vertex correction may be important 
since the Migdal theorem is not applicable. Although 
our results are consistent with the STS data, it is possi- 
ble that vertex correction may change a v hh , (q) . However 
the vertex correction to a v kk , (q) is beyond our scope in 
the present paper. 

We did not consider the el-ph matrix element between 
extended state and edge state. Although the matrix el- 
ement may be enhanced by the boundary deformation 
potential, it is naively expected that the matrix element 
is proportional to y '£,/ L and is negligible when L>(. In 
this case, the extended state and the edge state can be de- 
coupled. The geometry with d t ^S> L is referred to as the 
nano-graphite ribbon. For a ribbon, the off-site atomic 
deformation potential contributes to the scattering be- 
tween two edge states which are located at the different 
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edges since overlapping between the two edge states is 
not negligible. It is noted that Igami et al£^ showed 
that out-of-plane edge phonon modes appear depending 
on the effective mass of carbon atom at edge sites, and 
Tanaka et al*2» observed such modes at the armchair edge 
of nano-graphite ribbons on TiC(755) surface by high res- 
olution electron energy loss spectroscopy. Though b v As) 
and w„(q) used in this paper do not include the edge 
phonon mode, el-ph coupling for out-of-plane modes is 
negligible for the edge states. The el-ph interaction for 
nano-ribbon requires further studies on the el-ph inter- 
action. 

In summary, we formulated el-ph interaction for edge 
states and used it to calculate LDOS. Although our cal- 



culation does not include the Coulomb interaction, the 
result agrees with LDOS data£i&i when W « 0.1 eV. 
Our results should be compared with the future experi- 
ments of edge states in CNTs and graphene. 
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